Schwarz Type Solvers for hp-FEM Discretizations of Mixed Problems
نویسندگان
چکیده
The Stokes problem and linear elasticity problems can be viewed as a mixed variational formulation. These formulations are discretized by means of the hp-version of the finite element method. The system of linear algebraic equations is solved by the preconditioned Bramble-Pasciak conjugate gradient method. The development an efficient preconditioner requires three ingredients, a preconditioner related to the components of the velocity modes, a preconditioner for the Schur complementrelated to the components of the pressure modes and the discrezation by a stable finite element pair which satisfies the discrete infsup condition. The last condition is also important in order to obtain a stable discretization scheme. The preconditioner for the velocity modes is adapted from fast hp-FEM preconditioners for elliptic problems. Moreover, we will prove that the preconditioner for the Schur complement can be chosen as a diagonal matrix if the pressureis discretized by discontiuous finite elements.We will prove that the system of linear algebraic equationscan be solved in almost optimal complexity if the Qk − Pk−1,disc element is used. This yields to quasioptimal hp-FEM solvers for the Stokes problems and linear elasticity problems. The latter are robust with respect to the contraction ratio ν. The efficiency of the presented solver is shown in several numerical examples.
منابع مشابه
Inexact additive Schwarz solvers for hp-FEM discretizations in three dimensions
In this paper, a boundary value problem of second order in three space dimensions is discretized by means of the hp-version of the finite element method. The system of linear algebraic equations is solved by the preconditioned conjugate gradient method with an overlapping domain decomposition preconditioner with inexact subproblem solvers. In addition to a global solver for the low order functi...
متن کاملOn a parallel Maxwell eigensolver
Fast domain decomposition solver for internal problems of 3D hierarchical hp-FEM I. Anoufriev, V. Korneev Deriving fast domain decomposition (DD) preconditioners-solvers for hp discretizations of 2-nd order elliptic equations is still a challenge [1], although one pioneering solver of this type has been recently presented in [2] and generalized on adaptive hp discretizations in [3]. As it is we...
متن کاملPreconditioning High-Order Discontinuous Galerkin Discretizations of Elliptic Problems
In recent years, attention has been devoted to the development of efficient iterative solvers for the solution of the linear system of equations arising from the discontinuous Galerkin (DG) discretization of a range of model problems. In the framework of two level preconditioners, scalable non-overlapping Schwarz methods have been proposed and analyzed for the h–version of the DG method in the ...
متن کاملOn Higher Order Approximation for Nonlinear Variational Problems in Nonsmooth Mechanics
This paper is concerned with the hp-version of the finite element method (hp-FEM) to treat a variational inequality that models frictional contact in linear elastostatics. Such an approximation of higher order leads to a nonconforming discretization scheme. We employ Gauss-Lobatto quadrature for the approximation of the nonsmooth frictiontype functional and take the resulting quadrature error i...
متن کاملMultilevel Iterative Solution and Adaptive Mesh Reenement for Mixed Nite Element Discretizations
We consider the numerical solution of elliptic boundary value problems by mixed nite element discretizations on simplicial triangulations. Emphasis is on the eecient iterative solution of the discretized problems by multilevel techniques and on adaptive grid reenement. The iterative process relies on a preconditioned conjugate gradient iteration in a suitably chosen subspace with a multilevel p...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Comput. Meth. in Appl. Math.
دوره 12 شماره
صفحات -
تاریخ انتشار 2012